AMENDMENT 

In The Claims 

Please amend the claims as follows: 

1 . (Currently amended) An apparatus for computing multiple integral of a 
multidimensional integrand function A to be integrated with using a vector map mapping 
function f with unbounded support which converts an m (m>l)-dimensional vector having real 
number components into another m-dimensional vector having real number components, by 
which a multidimensional density function p for a limiting distribution resulting from repeatedly 
applying the map vector mapping function f to an m-dimensional vector u is analytically 
solvable, said apparatus comprising: 

a first storage unit which stores an m-dimensional vector x = (xj, X2, . . . ,x m ); 

a second storage unit which stores a scalar value w; 

a first computing unit which computes a vector x' = f(x) = (x'i, x'2, . . . , x' m ,) resulting 
from applying said vector map f to said vector x being stored in said first storage unit; 

a second computing unit which computes a scalar value w' = A(x)/p(x) based on said 
vector x being stored in said first storage unit and said scalar valu e w b e ing stor e d in said s e cond 
storag e unit ; 

an update unit which updates the value stored in said first storage unit by storing said 
vector x' computed by said first computing unit on said first storage unit, and updates the value 
in said second storage unit by adding said scalar value w' computed by said second computing 
unit to a value to be stored in said second storage unit; and 



2 



an output unit which computes, when the number of update times by said update unit 
becomes c (c>l), a scalar value s = w/(c+l) based on said scalar value w being stored in said 
second storage unit when the number of update times by said update unit becomes c (c>l) 5 and 
outputs said scalar value s as a result of the multiple integral. 

2. (Original) The apparatus according to claim 1, wherein said scalar value stored in 
said second storage unit first is a result from dividing a value resulting from applying said 
function A to said m-dimensional vector stored in said first storage unit first by a value resulting 
from applying said density function p to said m-dimensional vector stored in said first storage 
unit first. 

3. (Original) The apparatus according to claim 1, wherein said scalar value stored in 
said second storage unit first is 0; and 

said output unit computes a scalar value s' = w/c, and outputs said scalar value s 1 as the 
result of the multiple integral instead of said scalar value s. 

4. (Original) The apparatus according to claim 1 further comprising: 

a convergence rate obtainer which obtains convergence rate of scalar values sequentially 
output by said output unit while varying said number of update times c for each of plural vector 
maps gi, g2, . . . , gk (k>2) which are prepared as said vector map f; 
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a vector map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects a vector map gh (l<h<k) which shows fastest convergence rate; and 

an output controller which controls said output unit to output said scalar value with using 
said vector map gh as said vector map f and the number of update times c ! (c'>c) instead of said 
number of update times c. 

5. (Previously amended) The apparatus according to claim 1, wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 

u,f(u),f(f(u)),f(fl[f(u)))),... 

resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (ui 5 U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p(u)=ni=o m pi(ui); 
pi(Uj) ~ c_i|ui|' (1+a) for Uj -* — oo; 
pi(Uj) ~ c_i|ui|" (l+a) for Uj -» + oo; 
(a>0, l<i<m, c_j>0, c +i >0). 

6. (Previously amended) The apparatus according to claim 5, wherein said vector 
map f is defined as 

f(u)=(fi(Ui), f 2 (u s ), . . . , f m (Um)) 

by a function fj(t) = gj(di t)/d (dj>0) which is defined in l<i<m, and said map gi is 
defined by any one of the following maps q)j(l<j<8) and a natural number n\ (ni>2), as follows: 
g i ((p j (e))=( Pj (n i (e); 
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sgn(t) = 


= 1 for t >0; 


sgn(t) = 
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= -1 for t <(X 



7. (Currently Amended) The apparatus according to claim 5 further comprising: 
a convergence rate obtainer which defines said map f for each of plural positive numbers 

qi 5 q2 5 . . . , gk (^-2) prepared as an invariable a, and obtains convergence rates of the scalar 

values sequentially output by said output unit while varying said number of update times c; 
a positive number selector which refers the convergence rates obtained by said 

convergence rate obtainer, and selects an integer qh (l<h<k) which shows the fastest 

convergence rate; and 



an output controller which defines said map f with using said positive number qh said 
invariable a, and controls said output unit to output said scalar values with using the number of 
update times c' (c* >c) instead of said number of update times c. 

8. (Original) The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map gj with using plural ones of said 
maps <p>j, and obtains convergence rates of the scalar values sequentially output by said output 
unit while varying said number of update times c; 

a map selector which refers to the convergence rates obtained by said convergence rate 
obtainer, and selects one of said maps (pj which shows the fastest convergence rate; and 

an output controller which defines said map gj with using said map (pj selected by said 
map selector, and controls said output unit to output the scalar values with using the number of 
update times c f (cVc) instead of said number of update times c. 

9. (Original) The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map gi relating to each of plural natural 
numbers Pi, ?2, . . . , Pk (k>2) as said natural numbers n*, and obtains convergence rates of the 
scalar values sequentially output by said output unit while varying said number of update 
times c; 
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a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number Ph (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said natural number map gj with using said natural 
number Ph as said natural number nj, and controls said output unit to output the scalar values 
with using the number of update times c' (c'>c) instead of said number of update times c. 

10. (Original) The apparatus according to claim 1, wherein said output unit computes 
said scalar value s each time said update unit carries out update, compares said latest scalar value 
s with the former scalar value which is computed at former update, and outputs said latest scalar 
value s if a result of the comparison satisfies a predetermined condition for terminating the 
computation. 

1 1 . (Currently amended) A method for computing multiple integral of a 
multidimensional integrand function A to be integrated with using a vector map mapping 
function f with unbounded support which converts an m (m>l) -dimensional vector having real 
number components into another m-dimensional vector having real number components, by 
which a multidimensional density function p for a limiting distribution resulting from repeatedly 
applying the map vector mapping function f to an m-dimensional vector a is analytically 
solvable, said method using 

a first storage unit which stores an m-dimensional vector x = (xi, X2, . . . , x m ), and 



a second storage unit which stores a scalar value w, 
said method comprising the steps of: 

computing a vector x' = f(x) - (x'i, x f 2, . . . , 4) resulting from applying said vector map 
f to said vector x being stored in said first storage unit; 

computing a scalar value w ! = A(x)/p(x) based on said vector x being stored in said first 
storage unit and said scalar valu e w b e ing stor e d in said s e cond storag e unit ; 

updating the value stored in said first storage unit by storing said vector x computed by 
said first computing unit on said first storage unit, and updating the value in said second storage 
unit by adding said scalar value w' computed by said s e cond computing unit t o a value to be 
stored on said second storage unit; and 

computing, when the number of update times by said update unit becomes c (c>l), a 
scalar value s = w/(c+l) based on said scalar value w being stored in said second storage unit 
when the number of update times by said update unit becomes c (c>l), and outputting said scalar 
value s as a result of the multiple integral. 

12. (Original) The method according to claim 1 1, wherein said scalar value stored in 
said second storage unit first is a result from dividing a value resulting from applying said 
function A to said m-dimensional vector stored in said first storage unit first by a value resulting 
from applying said density function p to said m-dimensional vector stored in said first storage 
unit first. 
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13. (Original) The method according to claim 11, wherein said scalar value stored in 
said second storage unit first is 0, and 

said outputting step computes a scalar value s' = w/c, and outputs said scalar value s* as 
the result of the multiple integral instead of said scalar value s. 

14. (Currently Amended) The method according to claim 1 1 further comprising the 
steps of: 

obtaining convergence rate of scalar values sequentially output by said outputting step 
while varying said number of update times c for each of plural vector maps gi, g2 ? . . . , gk (k>2) 
which are prepared as said vector map f; and 

referring to the convergence rates obtained by said convergence rate obtainer, and 
selecting a vector map gh (l<h<k) which shows fastest convergence rate, and 

said outputting step outputs said scalar value with using said vector map gh as said vector 
map f and the number of update times c 1 (c'>c) instead of said number of update times c. 

15. (Previously amended) The method according to claim 1 1 , wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u))), . . . 

resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (ui, U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p(u)= iWpiCuD; 



Pi(uO~c_ i |u i r ,+,, > forui^-oo); 
Pi (uO ~ c_i |ui|- (,+a) forui^ + oo); 
(a>0, l<i<m, c-i>0, c+^0), 

16. (Previously amended) The method according to claim 15, wherein said vector 
map f is defined as 

f(u) = (f,(u 1 ),f 2 (u s ) 5 ...,f m (u m )) 

by a function fj(t) = gi(dj t)/dj (dj>0)) which is defined in l<i<m, and said map gi is 
defined by any one of the following maps q)j (l<j<8) and a natural number n h (n,> 2) as follows: 



gi <Pj(0)) = <Pj(ni0); 


9,(9) = 


-sgn(tan9) / |tan9| ,/a ; 


cp 2 (0) = 


-sgn(tan9) x |tan9| 1/a ; 


<P3(0) = 


-sgn(cos9) / |tan9| 1/a ; 


<P4(0) = 


-sgn(cos9) x |tan9| 1/a ; 


(p 5 (9) = 


sgn(cos9)/|tan9| 1/a ; 


<P6(9) = 


sgn(cos6)x|tan0| 1/a ; 


(p 7 (9) = 


sgn(sinO) / |tan9| 1/a ; 


<p 8 (9) = 


sgn(sin0)x|tan9| 1/a ; 


sgn(t) = 


= 1 for t >0; 


sgn(t) = 


= 0 for t =0; 


sgn(t) = 


= -1 for t <0, 



10 



17. (Original) The method according to claim 15 further comprising the step of: 
defining said map f for each of plural positive numbers qi, qz, . . . , qk (k>2) prepared as 

an invariable a, and obtaining convergence rates of the scalar values sequentially output by said 
output unit while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a positive number qh (l<h<k) 
which shows the fastest convergence rate; and 

defining said map f with using said positive number qh as said invariable a, and 
controlling output of said scalar values with using the number of update times c* (c*>c) instead of 
said number of update times c. 

18. (Original) The method according to claim 16 further comprising the steps of: 

defining said map gi with using plural ones of said maps cpj, and obtaining convergence 
rates of the scalar values sequentially output by said output step while varying said number of 
update times c; 

referring to the obtained convergence rates, and selecting one of said maps (pj which 
shows the fastest convergence rate; and 

defining said map g\ with using said selected map (pj selected, and controlling output of 
the scalar values with using the number of update times c* (c f >c) instead of said number of update 
times c. 
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19. (Original) The method according to claim 16 further comprising the steps of: 
defining said map gi relating to each of plural natural numbers Pi, P 2 , . . . , Pk (k>2) as 

said natural numbers n*, and obtaining convergence rates of the scalar values sequentially output 
by said outputting step while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a natural number Ph (l<h<k) 
which shows the fastest convergence rate; and 

defining said natural number map gi with using said natural number Ph, as said natural 
number n*, and controlling output of the scalar values with using the number of update times 
c' (c'>c) instead of said number of update times c. 

20. (Original) The method according to claim 11, wherein said outputting step 
computes said scalar value s each time said updating step caries out update, compares said latest 
scalar value s with the former scalar value which is computed at former update, and outputs said 
latest scalar value s if a result of the comparison satisfies a predetermined condition for 
terminating the computation. 

21 . (Previously amended) A computer readable recording medium storing a program 
for computing multiple integral of a multidimensional integrand function A to be integrated with 
using a vector map ma pping function f with unbounded support which converts an m (m>l) - 
dimensional vector having real number components into another m-dimensional vector having 
real number components, by which a multidimensional density function p for a limiting 
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distribution resulting from repeatedly applying the map vector mapping function f to an m- 
dimensional vector u is analytically solvable, said program comprising instructions for: 

storing an m-dimensional vector x = (xi, X2, . . . , x m ) in a first storage unit; 

storing a scalar value w in a second storage unit; 

computing a vector x 1 = f(x) = (x'i, x f 2, . . . x' m ), using a first computing unit, the vector x' 
resulting from applying said vector map f to said vector x being stored in said first storage unit; 

computing a scalar value w f = A(x)/p(x), using a first computing unit, the scalar value w' 
based on said vector x being stored in said first storage unit and said scalar value w b e ing stor e d 
in said s e cond storag e unit ; 

updating the value stored in said first storage unit, using an update unit, by storing said 
vector x' computed by said first computing unit on said first storage unit, and updates the value 
in said second storage unit by adding said scalar value w' computed by said second computing 
unit to a value to be stored on said second storage unit; and 

computing a scalar value s = w/(c+l), when the number of update times by said update 
unit becomes c (c>l), the scalar value s based on said scalar value w being stored in said second 
storage unit when the number of update times by said update unit becomes c (c>l), and outputs 
said scalar values as a result of the multiple integral. 

22. (Original) The recording medium according to claim 21, wherein said scalar value 
stored in said second storage unit first is a result from dividing a value resulting from applying 
said function A to said m-dimensional vector stored in said first storage unit first by a value 
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resulting from applying said density function p to said m-dimensional vector stored in said first 
storage unit first. 

23. (Original) The recording medium according to claim 21, wherein said scalar value 
stored in said second storage unit first is 0; and 

said output unit computes a scalar value s* = w/c, and outputs said scalar value s' as the 
result of the multiple integral instead of said scalar value s. 

24. (Original) The recording medium according to claim 21, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which obtains convergence rate of scalar values sequentially 
output by said output unit while varying said number of update times e for each of plural vector 
maps gi, g2, . . . , gk (k>2) which are prepared as said vector map t 

a vector map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects a vector map gh (l<h<k) which shows fastest convergence rate; and 

an output controller which controls said output unit to output said scalar value with using 
said vector map gh, as said vector map f and the number of update times c' (c*>c) instead of said 
number of update times c. 

25. (Previously amended) The recording medium according to claim 21, wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 
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u, f(u), f(f(u)), f(f(f(u))), 

resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (uj, U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p (u) = rU" 1 pi(ui); 
pKuO-c.ilUif 0 ^ forui^-oo); 
Pi (Ui) ~ c_i | Ui r C1+a) forui^ + oo); 
(a>0, l<i<m, c-j>0, c +i X)\ 



26. (Previously amended) The recording medium according to claim 25, wherein said 
vector map f is defined as: 

f(u) = (fi(ui), f 2 (u s ), . . . , f m (u m )) 

by a function fj(t) = gj(dj t)/dj (dj>0) which is defined in l^i^m, and said map gi is 
defined by any one of the following maps q)j (l<j<8) and a natural number nj (nj>2) as follows: 



gi (pj(0)) = (pj(ni8); 


9i(e) 


= -sgn(tan0)/|tan0| 1/a ; 


(p 2 (e) 


= -sgn(tan9) x |tane| 1/a ; 


93(6) 


= -sgn(cos9) / |tan0| 1/a ; 


(p 4 (9) 


= -sgn(cos0)x|tan6| 1/a ; 


<P5(6) 


= sgn(cos9)/|tan9| l/a ; 


96(6) 


= sgn(cos9)x|tan9| 1/a ; 


97(6) 


= sgn(sin9) / |tan0| 1/a ; 
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(P8(9) = sgn(sin9)x|tan0| l/a ; 
sgn(t) = 1 for t >0; 
sgn(t) = 0 for t =0; 
sgn(t) = -1 for t <0, 

27. (Original) The recording medium according to claim 25, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which defines said map f for each of plural positive numbers 
qi, q2, . . . , qk (k>2) prepared as an invariable a, and obtains convergence rates of the scalar 
values sequentially output by said output unit while varying said number of update times c; 

an integer selector which refers the convergence rates obtained by said convergence rate 
obtainer, and selects a positive number qh (l<h<k) which shows the fastest convergence rate; and 

an output controller which defines said map f with using said positive number qh, as said 
invariable a, and controls said output unit to output said scalar values with using the number of 
update times c 1 (c'>c) instead of said number of update times c. 

28. (Original) The recording medium according to claim 26, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which defines said map gj with using plural ones of said 
maps (pj and obtains convergence rates of the scalar values sequentially output by said output unit 
while varying said number of update times c; 
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a map selector which refers to the convergence rates obtained by said convergence rate 
obtainer, and selects one of said maps cpj which shows the fastest convergence rate; and 

an output controller which defines said map gj, with using said map (pj selected by said 
map selector, and controls said output unit to output the scalar values with using the number of 
update times c f (c'>c) instead of said number of update times c. 

29. (Original) The recording medium according to claim 26, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which defines said map gl relating to each of plural natural 
numbers Pi, P2, . . . , Pr (k>2) as said natural numbers nj, and obtains convergence rates of the 
scalar values sequentially output by said output unit while varying said number of update 
times c; 

a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number Ph (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said natural number map gj with using said natural 
number Ph as said natural number ni, and controls said output unit to output the scalar values 
with using the number of update times c' (c'>c) instead of said number of update times c. 

30. (Original) The recording medium according to claim 21, wherein said output unit 
computes said scalar value s each time said update unit caries out update, compares said latest 
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scalar value s with the former scalar value which is computed at former update, and outputs said 
latest scalar value s if a result of the comparison satisfies a predetermined condition for 
terminating the computation. 

3 1 . (Original) The recording medium according to claim 2 1 , wherein said recording 
medium is a compact disc, a floppy disk, a hard disk, a magneto-optical disk, a digital versatile 
disc, a magnetic tape, or a semiconductor memory. 
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